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Topological insulators |l|-[6|] is a new quantum 
phase of matter with exotic properties such as dis- 
sipationless transport and protection against An- 
derson localization. [7] These new states of quan- 
tum matter could be one of the missinglinks for 
the realization of quantum computing |8|, [9J and 
will probably result in new spintronic or magne- 
toelectric devices. Moreover, topological insula- 
tors will be a strong competitor with graphene 
in electronic application. Because of these poten- 
tial application the topological insulator research 
has literally exploded during the last year. Mo- 
tivated by the fact that up-to-date only few 3D 
systems are identified to belong to this new quan- 
tum phase |10l-ll8J we have used massive comput- 
ing in combination with data-mining to search for 
new strong topological insulators. In this letter 
we present a number of non-layered compounds 
that show band inversion at the T-point, a clear 
signal of a strong topological insulator. 

The prediction and realization of the 3D topological 
insulators 10l4l8| grew out of the work on t he q uantum 



Hall effect and quantum spin Hall effect in 2D [191-121] and 
the essence of these central articles [lll-l2l| is captured in 
a number letters describing topological insulators [l|, [3| — 

In school we have been taught that materials can be 
divided into insulators and conductors. Lately a new 
state of matter has emerged, the so-called strong topo- 
logical insulators in 3D. These materials can be viewed 
as insulators in the bulk but with metallic surface states, 
resembling the Dirac cones in graphene. These metallic 
surface states exhibit remarkable properties such as dis- 
sipation less transport and topological protection against 
perturbations and impurity scattering. Not surprisingly 
these new materials have spawned great attention and re- 
search in the science community because of the potential 
of interesting application. 

An insulator is characterized by its non-ability to con- 
duct, id est an external electric field does not cause cur- 
rent flow. Alternatively we think of an insulator as a 
system that lacks low lying excited states and it takes a 
finite amount of energy for an electronic excitation, i.e. 
we have an electronic band gap. It appears that the elec- 
tronic band gap is the one characteristic needed to char- 
acterize an insulator and this was true until the quantum 
Hall effect (QHE) came around. The quantum Hall state 
(QHS) shows up when, for example, a two-dimensional 
electron gas is subject to a magnetic field. The interior 



of the QHS has a band gap but contrary to a normal 
insulator the QHS has edge states that do not show a 
band gap, i.e. there is a current on the boundary of the 
sample. This very special charge flow is one directional 
thus making it protected against scattering [2[. Both 
the insulating- and QHS systems show bulk band gaps 
but the latter also have the surface states that are in- 
sensitive to scattering. So what is the difference between 
the two systems? It turns out that it is the topology of 
the occupied bands that does it, or more specifically, the 
topological class of the bulk band structure. 

Next came the quantum spin Hall effect (QSHE) 
20|, l21(. The quantum spin Hall state (QSHS) is char- 



acterized by yet another topological invariant (knot) and 
can exist without a magnetic field (a QHS required a 
magnetic field) thanks to the spin-orbit interaction. The 
bulk has a band gap where as the edge states do not. 
In graphene these edge states are spin-filtered thus mak- 
ing the different spins travel in opposite directions. The 
QSHE can give us edge states without a magnetic field 
in 2D thanks to the spin-orbit interaction. But what 
happens when we go 3D? In refs. Ill4l8| it was both 



predicted and experimentally verified that topological in- 
sulators do exist in 3D and it is the spin-orbit interaction 
that helps out. We can think of these edge states in 3D 
topological insulators as a higher dimensional version of 
the ID edge states in the QSHS. In 3D the simplest form 
of edge state can be viewed as a Dirac fermion metal 
with a linear excitation energy, c.f. the energy spectrum 
of mass less Dirac fermions (Dirac cone). Moreover, be- 
cause these surface states are topologically protected by 
time-reversal symmetry, Anderson localization does not 
occur even under strong disorder. A topological insu- 
lator is simply not allowed to lose the metallic surface 
state, id est become gaped or localize. This was recently 
experimentally verified for Bi x Sbi_ x [7[. 

When data-mining electronic structures in the search 
for new strong topological insulators the most straight 
forward method is to identify materials with a negative 
band gap [17|. To put it differently the electronic struc- 
ture should develop an anti-crossing feature at the T- 
point when the spin-orbit interaction is turned on as com- 
pared to no spin-orbit coupling. A typical band- inversion 
is shown in Figure [TJ Thus the data- mining algorithm is 
instructed to search for small gaped materials with an 
anti-crossing feature at the T-point. Note that the data- 
mine of electronic structures have been calculated with 
spin-orbit coupling. 

To show proof-of-concept the first three known sec- 
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FIG. 1: Band inversion because of strong spin-orbit coupling 
in small electronic band gap materials can be identified by 
searching for the anti-crossing feature around the T-point as 
the spin-orbit coupling is included in the calculation (right 
panel), c.f. without spin-orbit interaction (left panel). 



ond generation topological insulators Bi2Te3, Bi2Se3 and 
Sb 2 Te 3 



16Ml8| were investigated for band-inversion when 



spin-orbit interaction was turned on/off. Indeed all three 
materials develop the characteristic anti-crossing feature 
as spin-orbit coupling is included in the electronic struc- 
ture calculation. 
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FIG. 2: A schematic description of the search for new 
topological insulators. The key ingredient is the electronic 
structure project (ESP) that contain the electronic structures 
of tens of thousands of inorganic compounds. Data mining 
algorithms on the ESP has the potential to identify a variety 
of new functional materials. 



Before presenting results, the data-mining of electronic 
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Spgrp 


Struct. 


LDA 








type 


band gap 


[eV] 


Ca 3 PbO 


P m -3 m 


CaTiOs 


0.2 




Sr 3 PbO 


P m -3 m 


CaTiOs 


0.1 




Ba 3 PbO 


P m -3 m 


CaTiOs 


0.1 




Yb 3 PbO 


P m -3 m 


CaTiOs 


0.2 




Ca 3 SnO 


P m -3 m 


CaTiOs 


0.2 




Sr 3 SnO 


P m -3 m 


CaTiOs 


0.1 




Yb 3 SnO 


P m -3 m 


CaTiOs 


0.1 




GdPtSb 


F -4 3 m 


AlLiSi 


0.2 




Bi 2 SeTe 2 


R-3 mH 


Bi 2 Te 3 


0.3 




Bi 2 STe 2 


R-3 mH 


Bi 2 Te 3 


0.3 




PbTl 4 Te 3 


I 4/m c m 


In 5 Bi 3 


0.1 




BiTl 9 Te 6 


I 4/m c m 


In 5 Bi 3 


0.1 




BiTlTe 2 


R-3 mH 


NaCrS 2 


0.0 a 




SbTlTe 2 


R-3 mH 


NaCrS 2 


0.2 




Bi 2 TeI 


C 1 2/m 1 


Bi 2 TeI 


0.1 




GeSb 4 Te 7 


P -3 m 1 


AgBiSe 2 


0.2 




HgKSb 


P 63/m m < 


z KZnAs 


0.2 





a The material has small hole pockets. 

TABLE I: Results of the mining algorithm for second gen- 
eration non-trivial topological insulators. 17 compounds are 
identified as new potential topological insulators. ARPES 
measurements should confirm our findings. The electronic 
band structures, partial and total density of states as well as 
references to to the crystal structure determinations are given 
for all materials in the supplementary information. 



structures is briefly explained. The complete process 
is shown in Figure [2j The starting point is the inor- 
ganic crystal structure database (ICSD) which is a col- 
lection of some 130,000 experimental crystal structures 
obtained by X-ray and neutron diffraction. Using only 
the crystal structure as input we have calculated the elec- 
tronic structure for about 60,000 entries in the ICSD us- 
ing a full-potential linear muffin-tin orbital implementa- 
tion (FP-LMTO) [22] of density functional theory (DFT) 
within the local density approximation (LDA). The ob- 
tained library is maintained within the electronic struc- 
ture project (ESP) and has been made available on the 
web [23|. The data- mining process in the search for new 
strong topological insulators is rather straight forward. 
Non-layered small gap materials (<0.5eV) with an anti- 
crossing feature at the T-point are identified. Further- 
more it is verified that the anti-crossing feature disap- 
pears when the spin-orbit interaction is switched off. The 
data-mining identifies 17 compounds as potential strong 
topological insulators and these are presented in Table ID 
Note that because the ESP take crystal structures from 
the ICSD all compounds identified in this study exist and 
have been structurally determined. 

The first group of materials in Table U belong to 
the CaTiOs (perovskite) structure type but are inverse- 
perovskites. Perovskites are interesting materials be- 
cause these host a variety of materials properties such as 
superconductivity, spin dependent transport, ferroelec- 
tricity and colossal magnetoresistance. Present calcula- 



tions indicate that the strong topological insulator phase 
now can be added to this list. Through out the Bril- 
louin zone these materials have an in-direct LDA band 
gap ranging from 0.1 to 0.2 eV. The characteristic anti- 
crossing feature is present at the T-point and if the spin 
orbit coupling constant is set to zero the anti-crossing 
feature disappears. From the partial density of sates we 
conclude that these bands have mainly p-character. The 
electronic band structure for CasOSn is shown in Figure 

® 







FIG. 4: Bi2SeTe2 has the same structure type as the cele- 
brated topological insulator Bi2Te3- We see the same effect of 
the spin-orbit interaction as in the CasOSn case. The valence 
band at the T-point is pushed away and the anti-crossing 
feature show up when the spin-orbit interaction is included 
(black line, green line for no spin-orbit interaction). This 
compound could work under ambient conditions because of 
the 0.3 eV band gap. 



FIG. 3: CasOSn is an example of inverse perovskite that 
show band inversion. The material transforms from direct 
gap (no spin-orbit interaction included, green line) at the 
T-point to in-direct gap (around the same point in inverse 
space) when the spin-orbit is included (black line) in the 
electronic structure calculation, i.e. the upper most valence 
band is pushed way and the anti-crossing feature appears. 
The LDA band gap is 0.2 eV which qualifies the compound 
for operation under ambient conditions. 



Next in Table U is the half-Heusler compound GdPtSb 
with a AlLiSi structure type. The half Heusler family has 
been discussed in [10] but there with a distorted struc- 
ture. Here it is explicitly demonstrated for GdPtSb that 
with the experimental structure the material develops an 
anti-crossing feature when the spin-orbit term is included 
in the calculations and we speculate that this is also true 
for the other materials in [10]. GdPtSb has an indirect 
band gap of 0.2 eV. 

Bi2SeTe2 and Bi2STe2 both have the Bi2Te3 structure 
type. In fact these materials can be viewed as Bi2Te3 
but with every third Te replaced with Se and S, respec- 
tively. Bi2Te3 is a celebrated second generation topo- 
logical insulator [17| and both the Se and S variations 
have similar electronic structure to Bi2Te3 and show the 
anti-crossing feature at the T-point when spin-orbit in- 
teraction is removed. The indirect band gaps are 0.3 eV, 
respectively. The electronic band structure for Bi2STe2 
is shown in Figure HI For Bi2SeTe2 the electronic band 



gap closes without the spin-orbit term and the material 
appear to be a topological metal within the local den- 
sity approximation (LDA). However, LDA is known to 
underestimate the electronic band gap and using higher 
order corrections such as the GWA [24[ could very well 
open up a band gap thus making it a topological insula- 
tor. We speculate that there are many not yet discovered 
topological metals, i.e. a bulk metal but with topologi- 
cally protected surface states, and a data- mining search 
for these materials will be the subject for our next study. 

PbTUTe3 and BiTlgTe6 have a body-centered tetrag- 
onal Bravais lattice and belong to the InsBis structure 
type family. InsBi3 exhibit interesting electronic proper- 
ties such as superconductivity. The two compounds both 
show the anti-crossing feature as spin-orbit is included 
and we note the curvature of these bands is smaller com- 
pared to, for example, the Bi2Te3 family. The indirect 
band gap is 0.1 eV, respectively. 

SbTlTe2 and BiTlTe2 have the NaCrS2 structure type. 
The former material has an indirect band gap of 0.2 eV. 
The latter electronic structure is more complex. Using 
LDA the material show small electron and hole pocket 
but could very well turn out to have a small band gap 
when investigated using higher order theory such as 
GWA. 

The last three materials in Table [I] belong to different 
structure types and all have indirect band gaps but with 
more exotic electronic band structures. Bi2TeI develop 
the anti-crossing feature at two high symmetry points 
(r and V) in the Brillouin zone. GeSb 4 Te 7 and HgKSb 



show the band inversion at the A-point. The latter also 
display a small hole pocket away from the A-point. 

To conclude we have used massive computing and data 
mining to search for and identify new potential topolog- 
ical insulators. In Table U 17 new materials have been 
identified that all show an anti-crossing feature at the 
T-point (V, A for the last three) as the spin-orbit cou- 
pling is included. This feature is a clear signal for a 
topological insulator [17]. As pointed out in reference 
[10| it is crucial to identify new groups of materials (elec- 
tronic band structures) that are topologically non-trivial 
to maximize the likelihood for new physics as well as new 
devices @,[I1H. 



METHODS 

The more then 60,000 electronic structures used in 
the present study were calculated using a highly accu- 
rate full-potential linear muffin-tin orbital (FP-LMTO) 
[22| implementation of density functional theory (DFT) 
within the local density approximation (LDA) In essence 
DFT reduces the many-body problem of solving the elec- 
tronic structure of some 10 23 interacting electrons into to 
a one-body problem, thus making the calculation possi- 
ble. A more detailed description can be found in 22l. l23| . 



ACKNOWLEDGMENTS 

The Swedish Research Council (VR) and the Goran 
Gustafsson Stiftelse are acknowledged for financial sup- 
port. 



ADDITIONAL INFORMATION 



The authors declare no competing financial interests. 



* Mattias.Klintenberg@fysik.uu.se 
[1] Moore, J. E. The birth of topological insulators. Nature 

464, 194 (2010). 
[2] Hasan, M. Z. and Kane, C. L. Topological insulators. 

larXiv: 1002.38951 (2010): Colloquium article submitted to 

Rev. Mod. Phys. (2010). 
[3] Kane, C. L. An insulator with a twist. Nature Phys. 4, 

348 (2008). 
[4] Buttiker, M. Edge-state physics without magnetic fields. 

Science 325, 278 (2009). 
[5] Moore, J. Topological insulators; The next generation. 

Nature Phys. 5, 378 (2009). 
[6] Moore, J. An insulators metallic side. Nature 460, 1090 

(2009). 



[7] Roushan, P. et al. Topological surface states protected 
from backscattering by chiral spin texture. Nature 460, 
11061109 (2009). 

[8] Collins, G. P. Computing with quantum knots. Sci. Am. 
294, 57 (2003). 

[9] Bonderson, P. et al. A blueprint for a topologi- 
cally fault- tolerant quantum computer. Preprint at 
|http://arxiv.org/abs/1003.2856| (2010). 

101 Hsin, L. et al. Half-Heusler ternary compounds as new 

multifunctional experimental platforms for topological 

quantum phenomena. Nature Matr. 9, 546-549 (2010). 

I] Fu, L., Kane, C. and Mele, E. J. Topological insulators 

in three dimensions. Phys. Rev. Lett. 98, 106803 (2007). 

121 Moore, J. E. and Bals, L. Topological invariants of 
time-reversal- invariant band structures. Phys. Rev. B75, 
121306(R) (2007). 

131 R°y> R» Topological phases and the quantum spin Hall ef- 
fect in three dimensions. Phys. Rev. B79, 195322 (2009). 
11 Hsieh, D. et al. A topological Dirac insulator in a quan- 
tum spin Hall phase. Nature 452, 970 (2008). 

151 Hsieh, D. et al. Observation of unconventional quantum 
spin textures in topological insulators. Science 323, 919 
(2009). 

16] Xia, Y. et al. Observation of a large-gap topological- 
insulator class with a single Dirac cone on the surface. 
Nature Phys. 5, 398 (2009). 

171 Zhang, H. et al. Topological insulators in Bi2Se3, Bi2Te3 
and Sb2Te3 with a single Dirac cone on the surface. Na- 
ture Phys. 5, 438 (2009). 

181 Chen, Y. L. et al. Experimental realization of a three- 
dimensional topological insulator, BiTe. Science 325, 178 
(2009). 

191 Murakami, S. Nagaosa, N. and Zhang, S. C. Spin-Hall 
insulator. Phys. Rev. Lett. 93, 156804 (2004). 

201 Kane C. L. and Mele, E. J. Z2 topological order and the 
quantum spin Hall effect. Phys. Rev. Lett. 95, 146802 
(2005). 

211 Kane C. L. and Mele, E. J. Quantum spin Hall effect in 
graphene. Phys. Rev. Lett. 95, 226801 (2005). 
>1 Skriver, H. L. The LMTO Method: Muffin-Tin Orbitals 
and Electronic Structure. Springer, Berlin (1984). 

231 Klintenberg, M. http://gurks.fysik.uu.se/esp; Ortiz, 
C. Eriksson, O. and Klintenberg, M. Data mining and 
accelerated electronic structure theory as a tool in the 
search for new functional materials. Comput. Matr. Sci. 
44, 1042 (2009). 

241 Kotani, T. and Schilfgaarde, M. v Solid State Commun. 
121, 461 (2002). 

251 Fu, L. & Kane, C. L. Probing neutral Majorana fermion 
edge modes with charge transport. Phys. Rev. Lett. 102, 
216403 (2009). 

261 Qi, X-L. et al. Inducing a magnetic monopole with topo- 
logical surface states. Science 323, 11841187 (2009). 

27] Essin, A., Moore, J. E. & Vanderbilt, D. Magnetoelectric 
polarizability and axion electrodynamics in crystalline in- 
sulators. Phys. Rev. Lett. 102, 146805 (2009). 

281 Dzero, M. et al. Topological Kondo insulator. Phys. Rev. 
Lett. 104, 106408 (2010). 

29] Linder, J. et al. Unconventional superconductivity on 
a topological insulator. Phys. Rev. Lett. 104, 067001 
(2010). 

301 Ran, Y. et al. Composite fermions and quantum Hall 
stripes on the topological insulator surface. Preprint at 
|http://arxiv.org/abs/1003.09011 (2010). 



